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Observations of Large Scale Sidereal Anisotropy in 1 and 11 TeV cosmic rays from the MINOS 
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Abstract: The MINOS Near and Far Detectors are two large, functionally-identical, steel-scintillating sampling calorime- 
ters located at depths of 220 mwe and 2100 mwe respectively. The detectors observe the muon component of hadronic 
showers produced from cosmic ray interactions with nuclei in the earth's atmosphere. From the arrival direction of these 
muons, the anisotropy in arrival direction of the cosmic ray primaries can be determined. The MINOS Near and Far 
Detector have observed anisotropy on the order of 0. 1 % at 1 and 1 1 TeV respectively. The amplitude and phase of the 
first harmonic at 1 TeV are 8.2±1.7(stat.)x 10~ 4 and (8.9±12.1(stat.))°, and at 11 TeV are 3.8±0.5(stat.)x 10~ 4 and 
(27.2±7.2(stat.))°. 
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1 Introduction 

Cosmic ray primaries in the 1 TeV to 100 TeV energy 
regime are thought to originate within our galaxy H). 
Curvature in the galactic magnetic field (GMF) and multi- 
ple interactions with local scattering objects are expected 
to eliminate point-of-origin information and cause the 
observed arrival direction of cosmic ray primaries to be 
highly isotropic. However, large scale anisotropies at these 
energies have been observed by Milagro [2|, the Tibet 
Air Shower Array 0, ARGO-YBJ fW\, IceCube and 
Super-Kami okande 151 . 

Anisotropy can be introduced through multiple different 
mechanisms. Compton and Getting [6| predicted a dipole 
effect could arise due to the rotation of the solar system 
around the galactic center with the excess appearing in 
the direction of motion. Diffusion of cosmic rays out of 
the galactic-disk into the galactic-halo could introduce an 
anisotropy perpendicular to the disk (TJ. At TeV energies 
the heliosphere could produce a cosmic-ray excess in the 
direction of the heliotail; this excess should diminish with 
energy and is expected to be negligible above 10 TeV Q. 
Discrete sources have been shown capable of creating 
anisotropy at higher energies ifTOlfTTl . 

MINOS is a two detector long baseline neutrino beam 
oscillation experimentO- The 0.98 kton Near Detector, 
located at the Fermi National Accelerator Laboratory in 
Batavia, IL., is situated 94 m underground, below a flat 
overburden of 225 meters of water equivalent (mwe). The 



vertical muon threshold is 0.05 TeV. The raw atmospheric 
muon trigger rate is 27.29±0.01 Hz. The 5.4 kton Far 
Detector is located in the Soudan Underground Mine 
State Park in Soudan, MN., 720 m below the surface. 
The vertical overburden is 2100 mwe and the muon 
threshold is 0.73 TeV. The trigger rate at the Far Detector 
is 0.5364±0.0001 Hz. The MINOS detectors are two 
large, functionally-identical, steel-scintillating sampling 
calorimeters. The steel-scintillator planes are oriented ver- 
tically for optimal detection of horizontal beam neutrinos. 
As a consequence the detector has negligible acceptance 
for vertical muons increasing to near unity for horizontal 
muons. 

The MINOS detectors do not measure the arrival direction 
or energy of the cosmic ray primary. These quantities are 
inferred from the muons produced in the air shower. These 
muons are highly boosted objects and travel nearly parallel 
to the direction of the original primary. The mean surface 
energy for muons observed at the Near(Far) Detector is 
0.1 TeV(l.l TeV). The mean cosmic ray primary energy 
is ten times larger {8] than the mean muon surface en- 
ergy. The Near and Far Detectors therefore measure the 
anisotropy of 1 .0 TeV and 1 1 TeV primaries respectively. 

The expected level of cosmic ray anisotropy at TeV ener- 
gies is 0.1 %. Spurious sidereal anistropies can be created 
by the detectors' non-uniform sky exposure. The diurnal 
and seasonal variation in the muon flux due to variation 
in the atmospheric temperature can interfere to produce 
a anisotropy with a sidereal component. Both of these 
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effects must be removed to obtain a reliable measurement 
of the true cosmic ray anisotropy. The former has be 
mitigated by accurately measuring the detector live time; 
the latter by first measuring and then subtracting out the 
seasonal variation component. 



2 Event Selection 

The goal of this analysis is to identify large scale 
anisotropies in the arrival direction of atmospheric muons 
in the MINOS data set. The size of the data-set limits the 
search to anisotropies larger than 10°. These muons travel 
near parallel to the direction of the original cosmic ray 
primary; the large underground momentum of the muons 
limits the amount of multiple scattering in the overburden 
to be no more than 2 ° . The selection criteria have been 
chosen to obtain a detector pointing accuracy that, when 
combined with the scattering in the overburden and the 
initial scatter of the muon, gives an angular resolution 
on the direction of the primary cosmic-ray of 5°. Moon 
shadow searches performed using the MINOS Far Detector 
data have demonstrated that the detectors have a negligible 
pointing bias lfT2l . 

The same selections are applied to the data collected at the 
Near and Far Detectors. The candidate event must contain 
one and only one track. This track must be longer than 
2 m and possess a vertex which is no further than 50 cm 
from the detector edge. The x 2 /ndf of the track- fit must 
be less than 1, and the residual of the track to a straight 
line must be less than 0.04 m. Finally, it is required that 
the data be collected during a period of good detector 
run conditions. In summary the Far (Near) Detector had 
an effective live-time of 93.4% (76.0%), and a selection 
efficiency of 63.8 % (32.4 %). The Near Detector selection 
efficiency is much lower than that in the Far Detector as 
events in the down-stream end of the detector trigger the 
readout but are not reconstructed. After all selections we 
retain 67.7 x 10 6 events in the Far Detector and 0.989 x 10 9 
events in the Near Detector. 



3 Seasonal Variations 

When cosmic ray primaries interact with nuclei in the up- 
per atmosphere the subsequent showers contains kaons and 
pions. These mesons can either interact or decay producing 
the muons which we observe underground. The probability 
that a meson will decay is proportional to its energy and the 
density of the atmosphere through which it travels. Increas- 
ing the temperature of the atmosphere decreases the density 
reducing the probability that a meson will interact thereby 
increasing the observed muon rate underground. The rela- 
tionship between atmospheric temperature and muon rate 
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Figure 1 : The observed muon rate at the MINOS Near De- 
tector versus effective atmospheric temperature. The data 
is binned in 6hr intervals centered around the ECMWF 
data points. The errors on the muon rate are statistical. 
The temperature errors are set to ±0.5K. The red line is 
the best-fit of the data to equation [T] 



underground can be written as 
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where Ai? M is the deviation from the overall average muon 
rate. (Rn) is the muon intensity evaluated at a temperature 
(T c ff)- &t is the coefficient of correlation between rate 
and temperature. The temperature T e g- is an effective tem- 
perature of the atmosphere which weights the temperature 
of the atmosphere at different depths by the probability 
of a muon being produced at that elevation |fl3l . The 
temperature as a function of depth has been determined 
using the European center for Medium-Range Weather 
Forecasts(ECMWF) atmospheric model lfl"5l . The model 
provides atmospheric temperatures at 37 different pressure 
levels between 1 and lOOOhPa at four different times 
(0000 h, 0600 h, 1200h,1800h) throughout the day. 

The average muon rate and temperature at the MI- 
NOS Near (Far) Detectors is 8.7820±0.003 Hz 
(0.32566±0.00004 Hz) and 220.1±0.5K (222.2±0.5K) 
respectively. Figure [T] plots the observed muon rate at the 
Near Detectors versus the effective atmospheric temper- 
atures. To determine the value of olt a linear regression 
was performed on the data using ROOT'S MINUIT fitting 
package. The best fit values for ax are 0.421±0.004(stat.) 
for the Near Detector, and 0.886±0.009(stat.) for the Far 
Detector Q The value of ax is larger at the deeper Far 



1 . This measure of ar at the Far Detector is in agreement with 
our previous publication 1131 
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Detector as it samples from higher energy mesons which 
have an increased probability of interacting. 



4 Sidereal Anisotropy 

Detector acceptance effects dominate the number of events 
observed in declination making a true two-dimensional 
map, accurate to 0.1%, impossible. The rotation of the 
earth washes out these effects in right ascension. If no 
anisotropy exists, and provided we account for differing 
exposure as a function of detector right ascension^] we 
would expect the number of observed events to be flat as a 
function of right ascension. 

The rate-weighted number of muons as a function of right- 
ascension (a) and declination (S) are determined in each 1° 
detector right ascension bin k 
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where < R > is the average muon rate and R is the ex- 
pected muon rate at the current effective temperature T^. 
The two-dimensional sky-map is determined by summing 
over all detector right ascension bins and correcting for the 
different live times tj< in each bin. 
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where (i) is the average detector live time per detector 
right ascension bin. The level of anisotropy A in the right- 
ascension bin on for a given declination bin 8j is 

A(a- s . ) _ *M)-<m)> (4) 
A W> d >> <N(S J )> () 

where < N(6j) > is the average number of events 
per right ascension bin for that declination bin. The 
one-dimensional sky-map is obtained by integrating the 
data over all declination bins. 



A second order harmonic is typically used to describe the 
one-dimensional anisotropy 



deficit is observed in the Far detector data between 150° 
and 245°, the chance probability of this deficit occurring 
randomly is 0.001%. An 11.1 a deficit is observed in the 
Near Detector data between 155° and 225° with a chance 
probability much less than 0.0001%. The Near detector 
also sees a 7.8 a excess between 50° and 140° also with a 
chance probability much less than 0.0001%. 

The large statistics collected using the Near Detector are 
sufficient to allow an exploration of the two-dimensional 
anisotropy sky-map. Figure [3] plots the two-dimensional 
significance, defined as 



V< Wi) > 
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for 1 TeV CR primaries, where < N(Sj) > is the total 
number of events in the 5° smoothed data sample. A signif- 
icant deficit is observed Near the north galactic pole (192° 
ra, +27° dec), and an excess is observed around the helio- 
tail. The energy of our Near Detector data sample is equiv- 
alent to that collected by ARGO-YBJ experiment [ 16 1; and 
the two-dimensional sky-maps look similar. 



5 Conclusion 

The MINOS Near and Far Detector have been used to 
measure the anisotropy in the arrival direction of cos- 
mic rays at energies of 1 TeV and 1 1 TeV. The amplitude 
of the first-harmonic has been observed to decrease with 
increasing energy from 8.2±1.7(stat.)x 10~ 4 at 1 TeV to 
3.8±0.5(stat.)xl0~ 4 at 11 TeV. The two-dimensional sky- 
map at 1 TeV is consistent with previous publications. The 
MINOS experiment continues to take data and a more de- 
tailed analysis of the data will be forthcoming. 
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Figure 2: The one-dimensional projection of the Near and Far Detector sky-maps. The level of anisotropy in both data 
samples is 0.1%. The mean primary energy for the Far Detector is 1 1 TeV, and for the Near Detector is 1 TeV. 
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Figure 3: The two-dimensional significance sky-map for 1 TeV cosmic rays as determined using the MINOS Near Detec- 
tor. A 5 ° smoothing has been applied. 
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